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Abstract. A numerical investigation of two-laser photoassociation (PA) spectroscopy on spin-polarized 
metastable helium (He*) atoms is presented within the context of experimental observation of the least- 
bound energy level in the scattering potential and subsequent determination of the s-wave scattering 
length. Starting out from the model developed by Bohn and Julienne [Phys. Rev. A 60, (1999) 414], PA 
rate coefficients are obtained as a function of the parameters of the two lasers. The rate coefficients are 
used to simulate one- and two-laser PA spectra. The results demonstrate the feasibility of a spectroscopic 
determination of the binding energy of the least-bound level. The simulated spectra may be used as 
a guideline when designing such an experiment, whereas the model may also be employed for fitting 
experimentally obtained PA spectra. In addition, the prospects for substantial modification of the He* 
scattering length by means of optical Feshbach resonances are considered. Several experimental issues 
relating to the numerical investigation presented here are discussed. 

PACS. 34.20.Cf Interatomic potentials and forces - 32.80.Pj Optical cooling of atoms; trapping - 34.50.Gb 
Electronic excitation and ionization of molecules; intermediate molecular states (including lifetimes, state 
mixing, etc.) 



1 Introduction 

Within the context of ultracold quantum gases, Bose-Ein- 
stein condensation (BEC) of metastable helium has opened 
, up the way for experiments exploiting the 19.8 eV internal 
■ energy of the He* atoms in novel detection schemes ^ 
12131 • The accuracy of the results obtained in contempo- 
rary experiments, however, is compromised by the lack of 
an accurate determination of the value of the He* s-wave 
scat tering length |4T|o]|p]|7||g) ■ Another limitation in ultra- 
cold He* experiments is the absence of magnetic Feshbach 
resonances to control the interaction between the atoms, 
which is a consequence of the zero nuclear spin of 4 He. 

Standard techniques based on photoassociation (PA) 
spectroscopy have enabled accurate determination of the 
scattering length for most alkali elements. For the experi- 
mental determination of the scattering length two widely- 
used methods exist. The first one relies on determining the 
binding energy of the least-bound level in the scattering 
potential via a stimulated two-photon Raman transition, 
from which the scattering length can be inferred (see, for 
instance, Refs. This method is schematically de- 

picted in Fig.^ In the second method, Franck-Condon os- 
cillations observed in the linestrength of transitions from 
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the free-particle ground state to vibrational levels in an 
electronically excited molecular state are used to deter- 
mine the nodes in the scattering wavefunction, from which 
the scattering length may also be deduced 1 1 . 

For the case of two spin-polarized 2 3 S\ helium atoms, 
elastic scattering takes place in a 5 2J^ potential. However, 
the application of standard techniques to locate the least- 
bound state in this potential is not entirely straightfor- 
ward. The situation is complicated by the autoionization 
of metastable helium dimers, which occurs at short inter- 
nuclear range. Only dimers in a pure quintet spin config- 
uration at short range will experience suppression of au- 
toionization such that well-resolved vibrational resonances 
exist. This reduces a priori the number of useful pathways 
for a PA experiment. In addition, the short-lived singlet 
and triplet vibrational states effectively merge into a con- 
tinuum which may result in unwanted background signal 
during PA experiments. Indeed, long-lived vibrational lev- 
els have recently been observed in molecular states which 
possess quintet character at short internuclear range 
113) . Some of these states, which are located below the D-i 
(2 3 Si + 2 3 P2) asymptote, might be useful intermediate 
levels for a two-laser PA experiment aiming at the de- 
termination of the energy of the least-bound level in the 
scattering potential. 

An interesting alternative is offered by vibrational lev- 
els in the purely long-range 0+ potential well below Dq 
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Fig. 1. Schematic view of two-laser photoassociation, em- 
ployed here to convert a pair of free atoms (denoted as |/e,z)) 
into the least bound molecular state \g) supported by the scat- 
tering potential. 

(2 3 Si + 2 3 P ), observed by Leonard et al. [THT511T3] : 
see also Fig. |21 For such states the classical inner turn- 
ing point lies at long range, which leads to total inhi- 
bition of the short-range autoionization process. Despite 
their long-range location, the vibrational wavefunctions 
can be estimated to have sufficient Franck-Condon (FC) 
overlap with the least-bound wavefunction in the po- 
tential. The 5 £y~ potential has been calculated ab initio 
by Starck and Meyer ^B], from which a scattering length 
of 150 ao is inferred. Gadea et al. reconsidered the same 
potential and concluded that the least-bound vibrational 
level must be v = 14, with a binding energy of at least 
— 10 MHz jlj. This translates to a scattering length smaller 
than 340 ao - By comparison, the most accurate experimen- 
tal value for the He* scattering length to date (which was 
obtained by observing inelastic collisions at BEC thresh- 
old) is 195^20 a o U\- m agreement with the theoretical 
prediction. 

Control of the value of the scattering length may be 
achieved by means of an optical Feshbach resonance (OFR) , 
as first proposed by Fedichev et al. [TJ\, and recently ob- 
served by Theis et al. • It was pointed out that a laser 
coupling a colliding pair of ultracold atoms to a bound 
level in an excited molecular state may cause a significant 
modification of the scattering length. This scheme, how- 
ever, also implies loss of trapped atoms after spontaneous 
decay of the excited state. A suitable combination of ex- 
perimental parameters should be chosen to minimize these 
losses. 

Both two-laser PA and optical Feshbach resonances 
can be investigated numerically using a scattering model 
derived within the framework of generalized multichannel 
quantum defect theory, as done by Bohn and Julienne (19) . 
Their results, which are recapitulated in Section[21 include 
scattering matrix elements and a recipe for calculating 
two-laser PA rate coefficients. In Section l3~Tl two-laser PA 
spectra are simulated for the circumstances pertaining to 
the experiment described in Refs. (14115) . Next, in Sec- 
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tion l3.2l an account of the achievable modification of the 
He* scattering length using an OFR is given. In addition, 
experimental issues which may play a role in the choice 
of a suitable vibrational level for the OFR are discussed. 
Conclusions are presented in Section 0] 



2 Theoretical background 
2.1 Two-laser photoassociation 

Experimentally, PA spectra are obtained by monitoring 
the production rate of a certain product P resulting from 
the decay of the level being populated, while scanning 
the frequency of the PA laser. For instance, P could be 
the photon emitted during spontaneous decay. Plotting 
the rate coefficient associated with the production of P, 
Kp, versus PA laser frequency therefore simulates the PA 
spectrum. For a two-laser PA experiment involving a col- 
lection of atoms characterized by a Maxwell-Boltzmann 
energy distribution, Kp can be expressed as |2H] 

oo 

K P (T,Ui,U2,h,l2) = (2ttTiQ t )- 1 J2(21 + 1) 

1=0 

I* oo 

x / \S P {E,l,ui,U2,h,h)\ 2 e- E ' kBT dE, 
Jo 

(1) 

Here, the integration is over all initial scattering energies 
E. \Sp(E,l,uji,ui 2 , Ii, h)] 2 is the squared absolute value 
of the S- matrix element, which represents the probability 
of production of P as a function of the frequencies uj\ , u>2 
and intensities Ii,Js of lasers 1 and 2, respectively. The 
summation index I labels the partial wave contributing 
to the scattering process, and Qt = (/ifcsT/27r?i 2 ) 3 / 2 is 
the translational partition function with /i denoting the 
reduced mass of the system. For the two-laser case, there 
are three possibilities for the final state after the scat- 
tering process. These correspond to the transition of the 
initial free-atom pair (denoted as an energy-normalized 
state \Je.i)) to the excited state, |e), the bound molecu- 
lar ground state, \g), or back to |/) (see also Fig. For 
the case of PA of He*, the term "ground state" refers to 
the metastable (2 3 S + 2 3 S) 5 17+ state, and not the ac- 
tual ground state of the helium dimer. The probability for 
each of these transitions to take place is determined by the 
squared absolute value of the corresponding S'-matrix ele- 
ment, the explicit form of which is given in Eqs. (4.8-4.11) 
of Ref. JU] . The S'-matrix elements themselves depend on 
the rate 7 S of stimulated free-to-bound transitions induced 
by laser 1, the bound-bound Rabi coupling by laser 2, fi, 
the decay widths r e and r g of the excited and ground 
molecular bound states, respectively, and the detunings 
Ai and A^ of lasers 1 and 2. The detunings are defined 
as Ai — E e /h — u>i and A% = E g /h — (a>i — u^), where E e 
and Eg are the energies (with respect to the 2 3 S + 2 3 S 
asymptote) of the excited-state and ground-state vibra- 
tional levels, respectively. This definition of the detuning 
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implies that a negative value for Ai corresponds to a blue 
detuning. The stimulated rate 7 S is obtained from Fermi's 
golden rule: 



l s =(2-K/%)\{f E j\V raA \e)Y 



(2) 



Vrad = d 



h_ 

e c 



1/2 



(3) 



Here, d is the atomic dipole moment operator, and the po- 
larization of the light is ignored for the moment. It should 
be mentioned that the use of energy-normalized states au- 
tomatically includes the correct Wigner-law threshold be- 
haviour, which means that 7 S vanishes as E^ 2l+1 ^ 2 for 
E -> 0. 

Vrad also determines the Rabi coupling between the 
ground and excited bound states, and the Rabi frequency, 
f2, is given by 

= {e\V tad \g)/%. (4) 



2.2 Optical Feshbach resonance 

To create an optical Feshbach resonance, only one PA laser 
is required. The S-matrix elements for two-laser PA may 
also be used to describe PA by a single laser by setting 
the intensity of laser 2 to zero. Then, from the scattering 
matrix element Soo for elastic scattering (see Eqs. (4.9- 
4.11) in Ref. ^5]) a complex- valued phase shift, 5, may be 
obtained using the relation (23 



Soo = CX P (2i5) • 



(5) 



For low scattering energies, this complex phase shift trans- 
lates to a complex- valued scattering length, a, according 
to ( 

a = A — iB = — lim — tan S. (6) 



fe^o k 



Here, k 



2fj,E/h 2 denotes the initial wave number. 



Using the explicit expression for Soo given in Ref. |19| . it 
can be shown that 



A = a nat + 



1 ls [E - hA 1 - E{h)] 
2k [E-hA 1 -E{I 1 )] 2 + (r e /2) 



(7) 



where a nat is the natural value of the scattering length, 
and 

i 7s r e 

B= 4k[E-HA 1 - E{h)] 2 + (T e /2) 2 ■ (8) 

Here, E(Ii) is an energy- and intensity-dependent fre- 
quency shift ^5], which will be ignored below as it does not 
fundamentally change the results. The imaginary part of a 
accounts for losses due to spontaneous emission occurring 
in the excited state. By tuning the intensity and/or the 
detuning from the resonance, a may be varied. It should be 
noted that within this picture, mixing due to state dress- 
ing of the initial free state and the excited bound state 
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where the radiative coupling, y ra d , by laser 1 is given by 3 
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Fig. 2. PA of spin-polarized He* atoms. Laser 1 (dashed ar- 
row) can couple two free atoms (which are initially in the Sg 
scattering potential) to an excited bound state, for instance in 
the potential correlating to the D2 asymptote, thereby creat- 
ing an OFR. The dashed short-range part of the potential well 
indicates that it is not accurately known. Laser 1 may also be 
tuned to v = vibrational state in the purely long range 0i 
potential below the Do asymptote (solid arrow). Laser 2 may 
then resonantly couple this state to the least-bound (v — 14) 
vibrational state in the scattering potential. The differ- 

ence frequency of the two lasers is directly related to the bind- 
ing energy of the v = 14 state. 



is ignored. This has been predicted to give rise to addi- 
tional frequency shifts and even extra resonances for PA 
at relatively high intensity |22II23| . 

The real and imaginary parts of a can be converted into 
rate coefficients, K e \ and -Kinelj for elastic and inelastic 
collisions, respectively |21) : 



K el = k(A 2 

A 1 



B 



incl 



8tt?i 



B. 



(9) 



Here, i^inei relates to scattering events in which two atoms 
are lost (in Section HOI the various loss processes will be 
described) . 



3 Numerical results 

3.1 Two-laser photoassociation spectra 

As has been pointed out in Ref. ^3], the v = 0, J = 1 
state in the 0+ potential well below Do (Fig. [2J is ex- 
pected to have sufficient FC overlap with the least-bound 
state in the 5 17+ scattering potential. The value of the 
squared FC-overlap integral is of the order of 10 _1 (which 
is more than one order of magnitude larger than for any 
of the other five bound states observed in the 0+ well), 
for which case the efficiency of the second step in the two- 
laser PA process should be reasonable. The fact that the 
v — 0, J = 1 line could be observed with a good signal- 
to-noise (S/N) ratio indicates that also the the first step 
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can be made at a sufficiently high rate For the sim- 
ulation of two-laser PA, the v = 0, J = 1 state will be 
used as the intermediate state, and conditions similar to 
those in the experiments described in Refs. |14ll5j will be 
assumed. These imply PA on a magnetostatically trapped 
cloud consisting of 10 5 -10 6 spin-polarized He* atoms, at 
densities 10 12 -10 13 cm -3 , and temperatures 10-100 /iK. 
It was checked that for such low temperatures, all I > 2 
contributions to the summation in Eq. are negligibly 
small as compared to to the s-wave contribution. Since 
odd partial waves do not contribute to the scattering of 
identical bosonic He* atoms, only the / = term is re- 
tained. This is consistent with the fact that no vibrational 
levels in 0+ with total rotational quantum number J > 1 
were observed experimentally fT4"HT5] . 

For the evaluation of the matrix elements in Eqs. J5J 
and (0J, the wavefunctions of the involved states must be 
known. For the 5 ^ wavefunctions, these are obtained 
by numerically solving the radial Schrodinger equation 
using the potential data from Ref. [To]. Of course, the 
wavefunctions depend on the actual shape of the poten- 
tial, which is not accurately known. For the simulation of 
two-laser PA spectra the above mentioned potential data 
are used. In case that the actual potential turns out to 
deviate from the one used here, the difference between 
the two potentials will predominantly affect the radia- 
tive coupling (Eq. ©) via the FC overlap of the radial 
wavefunctions (see also the Appendix). However, for most 
potential curves lying within the range of the theoretical 
prediction, the differences due to this effect will not be 
too large. Also the 0+ wavefunction is found by numerical 
integration of the Schrodinger equation, using the Born- 
Oppcnheimer potential obtained from diagonalization of 
the long-range interaction between a 2 3 S and a 2 3 P atom. 
Here, the long-range interaction includes the electrostatic 
resonant dipole-dipole interaction (which scales with the 
internuclear distance R as R~ 3 ), the Van der Waals in- 
teraction (scaling as R~ 6 ), and the atomic relativistic in- 
teractions (spin-orbit, spin-spin etc.) |15| . Further details 
of the evaluation of the matrix elements are postponed to 
Appendix B. 

Having acquired all necessary information, Eq. JQ) can 
be employed to calculate the rate coefficient Kf for ob- 
servation of fluorescence decay of the excited state. In 
Fig. Ela), K f is plotted versus laser-frequency detuning 
for a single-laser PA experiment, and for various temper- 
atures. For the 25 /LtK case, the line profile has a width 
close to the spontaneous decay width of ~ 3 MHz, which 
is about twice the width of the atomic 2 3 Si — > 2 3 P2 
line |15| . For higher temperatures the line becomes in- 
creasingly more shifted, broadened and skewed. It also 
follows that for intensities well below 10 4 7 sa t (with J sat = 
1.67 W/m 2 the saturation intensity of the atomic 2 3 Si — > 
2 3 P2 transition), broadening due to saturation does not 
play a significant role. The bandwidth of the PA laser is 
neglected, which is justified for a laser with a few-hundred 
kHz bandwidth. 

For the two-laser case depicted in Fig. |3Jb), laser 1 
is scanned across the Autler-Townes doublet induced by 
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Fig. 3. Rate coefficients for observation of spontaneous emis- 
sion from the excited state plotted as a function of Ai, for 
various temperatures, (a) Single-laser lineshapes with Ii — 
100 7 sat . (b) Two- laser lineshapes, with A 2 = 0, I\ = 100 7 sa t 
and I2 = 400 J S at • The spectra in (b) correspond to the case of 
frustrated PA. 



laser 2. Again, at higher temperatures the spectrum be- 
comes more shifted, and the Autler-Townes doublet be- 
comes less resolved. In Fig. 0{a), an example of stimu- 
lated two-photon Raman PA is depicted. Laser 2 is now 
detuned several lincwidths below resonance, and scanning 
laser 1 from blue to red shows the " unperturbed" line and 
the narrow Raman line, respectively. The width of the lat- 
ter is determined by the autoionization width of the 5 ^^" 
state (see below), and by thermal broadening. An alter- 
native approach for detecting the two-photon transition is 
shown in Fig.^Jb). Here, laser 1 is kept at the maximum 
of the single-laser line (i.e. A\ « 0), whereas laser 2 is 
scanned. As long as laser 2 is not resonant, a background 
signal due to PA will be present. When laser 2 becomes 
resonant, the single-laser line will split up into the Autler- 
Townes doublet, causing a reduction of PA signal. This 
"frustrated PA" scheme may be convenient for an initial 
search, as the decrease in PA signal can be made 100%. 

At this point, the connection with the experiment de- 
scribed in Refs. |14II15| demands some attention. In this 
experiment, the PA spectra were obtained in a magnetic 
trap using a novel calorimetric detection scheme, which 
is briefly outlined below. Spontaneous decay of the ex- 
cited state will produce fast atom pairs with energies up 
to ~ 30 mK, with some of the pairs being in trapped 
magnetic substates. If the temperature of such atom pairs 
is below the 10 mK trap depth, their kinetic energy will 
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Fig. 4. Rate coefficients for observation of fluorescence 
decay (Kf, solid curve) and for observation of autoionization 
of 2 3 S + 2 3 S molecules formed by stimulated two- photon 
transitions (Kai, dotted curve), (a) Two- laser lineshapes with 
h = 100 /sat, h = 400 / sa t, and A 2 = 8 MHz. (b) Frustrated 
PA with Ii — I2 = 100 /sat- Here, Z\2 is scanned while zii ~ 
(corresponding to the peak of the single-laser line, whose height 
is indicated by the dashed horizontal line) . The temperature is 
25 for all plots 



be deposited in the trapped cloud by rethermalizing colli- 
sions. For a sufficiently cold cloud, this leads to substantial 
heating due to resonant PA, while the decrease in atom 
number due to the escaping atom pairs may also be de- 
tected. Since the heating is related to spontaneous decay 
of the excited state, the spectra plotted in Figs.[3]may also 
resemble those measured using the calorimetric method. 
However, although it can be estimated that the above 
mechanism nearly explains all the heating in the one-laser 
PA experiments [23], other mechanisms should be con- 
sidered too. Spontaneous decay of the excited state may 
also result in a bound S'g state, which will eventually 
decay by autoionization. Due to its small binding energy, 
the radial wavefunction of the v = 14 vibrational state 
extends far enough to have reasonable FC overlap with 
the purely long range vibrational states in J , which have 
classical inner turning points near 150 ao. By contrast, the 
outer turning points of the vibrational states in 5 with 
v < 13 lie at internuclear distances smaller than 30 ao |16| . 
Therefore, the v = 14 state is most likely to be populated 
via spontaneous decay. Combining a semiclassical estimate 
of the v = 14 vibration frequency with the Penning ion- 



ization suppression factor given in Ref. |25| . one arrives 
at an approximate value for the inverse autoionization- 
limitcd lifetime r g = (5 /is) -1 . The autoionization will 
produce fast He + or HeJ ions, which may cause heating 
via collisions with trapped He* atoms similar to the heat- 
ing observed in a He* BEC 0. This process may become 
important in the case of stimulated production of 
molecules. To illustrate this, Figs. 0] display the rate coef- 
ficients Kai for decay by autoionization. Heating due to 
Kai should be accounted for by multiplying i^Ai with a 
certain constant and adding it to Kf to simulate the calori- 
metric spectrum. Likewise, for the frustrated PA spectrum 
depicted in Fig.0fb), the heating due to 5 E+ bound state 
formation will shallow the dip near A2 = 0. However, as 
Kai for this case is more than two orders of magnitude 
smaller than Kf , this effect might be relatively small. The 
lineshapes for Kai in Figs. 0] also suggest the interesting 
option of ion detection for observing the two-laser transi- 
tion. Finally, heating due to off-resonant atomic scatter- 
ing should be taken into account. For intensities of a few 
hundred I sa t, the off-resonant atomic scattering rate can 
be calculated to be of the same order of magnitude as 
the PA rate for a cloud at 10 13 cm -3 density. Although 
the photon recoil associated with each atomic scattering 
event contributes to the heating by about twice the recoil 
energy (which is small compared to the heating per molec- 
ular decay event), the scattering may leave the atom in a 
different magnetic substate for which Penning ionization 
is not suppressed. A subsequent Penning-ionizing collision 
may cause significant heating of the cloud. On the other 
hand, heating due to off-resonant scattering varies slowly 
with detuning compared to heating due to PA. In a PA ex- 
periment, If may be chosen so as to obtain an sufficiently 
large rate coefficient Kf with a minimum of background 
signal, thus leading to optimum S/N ratio. 

Apparently, the PA signal obtained from heating or 
loss is related to the strength of the PA process only in 
an indirect way. This makes the calorimetric detection 
method less suitable for observing FC oscillations as in 
Ref. PH. 

The model used here to simulate two-laser PA spec- 
tra may also be employed for fitting experimental data. 
With the temperature known from experiment, the inten- 
sity of laser 2 and the binding energy of the 5 £y~, v = 14 
level may be used as fit parameters. Moreover, the He* 
2 3 S + 2 3 P system has good theoretical access: radiative 
coupling can be computed using accurate theoretical data, 
and hyperfine interaction and higher-order partial wave 
contributions are absent. The two-laser PA spectroscopy 
proposed in this work may thus serve as an absolute test 
case for existing theoretical models. 



3.2 Optical Feshbach resonance 

An OFR requires relatively intense laser fields, which may 
cause heating or loss of atoms due to off-resonant atomic 
scattering. Fedichev et al. pointed out that choosing a vi- 
brational level with large binding energy reduces the losses 
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due to atomic scattering, despite the fact that higher in- 
tensities will be needed to achieve the same modification 
of a 17] . The purely long range states below Dq discussed 
above have comparatively small binding energies with re- 
spect to the D Q asymptote, and the laser intensity re- 
quired for an OFR can be estimated to cause strong off- 
resonant atomic scattering. The subsequent heating and 
losses would severely limit the lifetime of a cold cloud of 
atoms, and preclude the use of the purely long range states 
for an OFR. 

As an alternative, the present calculation assumes the 
existence of a vibrational level with a binding energy of 
-250 GHz, located in a (2 3 S + 2 3 P) 5 27+ potential be- 
low D 2 - The 5 27+ potential couples at long-range to the 
Hund's case (c) states 1„ and 2 U , for which is the inter- 
action is well known |26U15j . The short-range part of this 
potential has never been considered by experiment or the- 
ory. Nevertheless, it is to be expected that autoionization 
of this state is strongly suppressed, similar to the sup- 
pression of Penning ionization during collisions between 
spin-polarized He* atoms J3| ■ The long-range part of the 
vibrational wavefunction is obtained by inward numerical 
integration of the radial Schrddinger equation, assuming 
a binding energy of —250 GHz and taking only the long- 
range part of the interaction into account. For simplicity, 
the long-range part is represented by a —C3/R 3 — Cq/R 6 
potential (with the dispersion coefficients from Ref. 26 ), 
which has a negligible effect on the wavefunction of the 
state bound by —250 GHz. The inner part of the excited- 
state wavefunction thus remains unknown. However, this 
part of the wavefunction oscillates rapidly with R, and it 
can be estimated to contribute only marginally to the FC 
overlap integral involved in the matrix element in Eq. (J5J) • 
Now, the effect of PA light (nearly resonant with the vibra- 
tional level) on the scattering length and the inelastic loss 
rate coefficient can be studied using Eqs. 10, ©, and ©. 
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Fig. 5. Optical Feshbach resonance using a PA intensity of 
5 x 10 5 W/m 2 . Solid curve: real part of the scattering length as 
a function of detuning away from a vibrational resonance with 
binding energy —250 GHz. The natural value of the scattering 
length, a nat , is indicated by the dotted horizontal line and is 
assumed to be 200 ao. Dashed curve: rate coefficient for loss of 
atoms due to spontaneous decay of the excited state. 
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Fig. 6. Parametric plot of Kin e \ versus A, for several PA laser 
intensities and with Ai being varied as a parameter. Here, a na t 
is assumed to be 200 ao. 

Fig- El shows the familiar dispersive behaviour of A versus 
detuning (A may even become negative) , and the resonant 
character of Ki ne \. Inelastic losses stem from decay of the 
excited vibrational level by spontaneous emission, which 
may result either in fast and/or untrapped atoms or in 
autoionizing 2 3 S + 2 3 S molecules, or from autoionization 
of the excited state at short internuclear range. 

A more convenient graphical representation of the re- 
lation between A and K lne \ is shown in Fig. [S] Here, A 
and i"Li nc i are plotted for various PA intensities, with the 
detuning being varied as a parameter. It is clear that a 
larger intensity allows for larger modification of the scat- 
tering length at the same level of losses. For the particular 
case in which A is made to vanish, the losses are seen to 
decrease with increasing PA intensity (here it should be 
mentioned that for increasing PA intensity, A becomes 
zero for increasing values of Ai). Another crucial role is 
played by the natural value of the scattering length, a na t . 
This is illustrated in Fig. da), where the maximum achiev- 
able value for A is plotted as function of PA intensity for 
several values of a na t lying within the uncertainty of the 
theoretical prediction [27J . It follows that for a larger value 
of a n at the modification will be larger too. Similarly, the 
parameters determining a vanishing value for A depend 
strongly on a na t as well. As a consequence, the inelastic 
loss rate coefficient for the case A = decreases with in- 
creasing values of a nat , as can be seen in Fig. Efb). This 
strong dependence on a nat emphasizes the importance of 
an accurate (spectroscopic) determination of the scatter- 
ing length. 

For any of the cases discussed above, the timescale 
on which the temperature and number of atoms in the 
cold cloud will be seriously affected is small. For a cloud 
of He* atoms at a density of 10 13 cm' 3 , PA losses set 
a timescale of 10 _5 -10 -3 s, whereas off-resonant atomic 
scattering typically limits the lifetime to several millisec- 
onds or longer. These figures confirm the results of an ear- 
lier study of prospects for optical Feshbach resonances in 
He* |28j . In fact, the elastic rate never exceeds the inelas- 
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Fig. 7. (a) Maximum achievable value for A versus PA inten- 
sity I\. The various curves correspond to different values for 
dnat- (b) Inelastic loss rate coefficient for the case of vanishing 
A versus PA intensity, for several values of a nat . 



tic rate by far, and it even is smaller in most of the cases 
discussed above. This means that the present scheme for 
optical control of the scattering length can be used only 
for applications in which short, impulse- like changes of the 
scattering length are required. For studies of a He* BEC 
in an optical lattice, which offers the unique feature of de- 
tection of ions produced in Penning collisions, tunability 
of the scattering length may be interesting [2H|. For in- 
stance, for the case of two He* atoms occupying the same 
lattice site, the PA rate (and thus the modification of a) 
may be enhanced [3U] ■ This would reduce the losses associ- 
ated with the conditions which set a to a particular value. 
Another point of further study may be the behaviour of 
the Penning-collision rate as a function of a. 

When designing an experiment to optically control the 
scattering length, the binding energy of a suitable vibra- 
tional level in the excited-state potential should be known. 
In the case of alkali elements, the location of suitable 
bound states can often be obtained from existing spec- 
troscopic data on the excited-state potential. In contrast 
to the case of alkali atoms, helium dimers do not natu- 
rally exist, which explains the poor amount of spectro- 
scopic data on the 2 3 S + 2 3 P states. Theoretical data 
are available only for the (2 3 S + 2 3 P) 3 17+ state; data 
of singlet and quintet molecular potential energy curves 
have hitherto not been published. To avoid an extensive 
spectroscopic search for lines near the desired binding en- 
ergy, the accumulated-phase method might be employed 
to assist in predicting their locations |31II12| . For an accu- 
rate prediction of vibrational level positions, however, the 
availability of ab initio data is probably crucial. 

When deciding on a suitable vibrational level to use 
for an optical Feshbach resonance, a point of concern is 
whether nearby molecular vibrational resonances in sin- 
glet or triplet potentials can be excited by the PA light. 
Such resonances, which will be broadened to almost the 



vibrational level spacing because of the large (near unity) 
probability of autoionization at short range |32| , introduce 
additional losses. However, for the Condon radius corre- 
sponding to excitation at —250 GHz, all nearby molecu- 
lar potentials are nearly pure Hund's case (a) states, as 
the relativistic atomic interactions that mix them at long 
range are small compared to the resonant dipole-dipole 
interaction. Starting from a sample of spin-polarized He* 
atoms, this means that only quintet ungerade states can 
be accessed, and no PA to molecular states coupling to 
strongly autoionizing singlet and triplet configurations is 
possible. 



4 Conclusion 

Two-laser PA spectroscopy of ultracold spin-polarized He* 
atoms, aiming at the detection of the least-bound vibra- 
tional level in the 5 7>7 + scattering potential, has been in- 
vestigated numerically. The binding energy of this state 
may be used to infer an accurate value for the He* s- 
wave scattering length. The scheme, which assumes the 
purely long range v = vibrational level in the 0+ poten- 
tial well below Dq as the intermediate level, may result in 
a detectable two-photon stimulated Raman transition for 
feasible experimental parameters. Different approaches to 
observe the two-laser transition are discussed within the 
context of a calorimetric detection method |14U15j . 

In addition, optical control of the He* scattering length 
using optical Feshbach resonances has been explored nu- 
merically. Modification of the scattering length to very 
large (> 10 3 ao) values or to vanishing or even negative 
values is anticipated for experimentally realizable condi- 
tions, although the magnitude of the optical Feshbach res- 
onance depends strongly on the natural value of the scat- 
tering length. Also, this method to control the scatter- 
ing length may be applied only for a very short (10~ 5 - 
10 -3 s) duration, to avoid substantial loss and/or heating 
of trapped atoms. This will undoubtedly complicate a pos- 
sible experimental implementation of an OFR. Notwith- 
standing all the experimental disadvantages as compared 
to magnetic Feshbach resonances, optical Feshbach reso- 
nances are presumably the only handle one has to modify 
the value of the He* scattering length. Finally, several ex- 
perimental issues concerning the search for and choice of 
a suitable excited-state vibrational level have been dis- 
cussed. 
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6 Appendix: matrix elements of the radiative 
coupling 

The matrix element in Eq. (J2J involves the energy-norma- 
lized state (the label I is dropped since only the I = 
contribution is taken into account), which can be written 
as the product of a radial wavefunction, a rotationless an- 
gular part, and the rotational part of the molecular state: 



also depends on R via the expansion coefficients Ci(R) 
in Eq. (|12J) . The FC integral is therefore evaluated over 
the full internuclear range of overlap, including the radial 
wavefunctions and expansion coefficients Ci(R), with each 
coefficient weighted by the value of the corresponding an- 
gular part of the matrix element. The latter is evaluated 
using the direction cosine matrix E3 , 34 . It should also be 



1 S,e{R) 
R 



5 E+)\Jf2M). 



(10) 



Here, fl is the projection of J on the internuclear axis, 
whereas M is the projection of J on the lab-fixed z-axis, 
defined by the field of the magnetic trap. For spin-polarized 
He* atoms, M = 2; hence J = 2 (for I = 0). Q de- 
pends on the alignment of the internuclear axis with re- 
spect to the z-axis. The energy-normalized radial wave- 
function, Uf t E(R), is obtained by solving the Schrodinger 
equation for two free atoms in the scattering potential 
with an asymptotic energy E, with the long-range bound- 
ary condition [B?] 



u fM R ) 



2fi 
Trh 2 k 



sm[kR + rj(k)], R -> oo. (11) 



Here k is the initial wave number, and r)(k) denotes an 
energy-dependent phase shift. The left-hand side of Eq. I|ll|) 
corresponds to an energy-normalized, free-particle state. 
The excited bound state can be written as 

|e) = Uv= i R) \Q+)\J'n'M') 



R 

u v=0 (R) 



{ Cl ( J R)rr+)+c 2 (i?)| 6 iT„) 



R 

|n.c.)}|J'n'M'), 



(12) 



where Ci(R) are the expansion coefficients of |0+) expressed 
in ungerade Hund's case (a) states. |n.c.) is the sum of all 
singlet and triplet states in the expansion |0+), which are 
not coupled to the 5 17+ lower state by a dipole transition. 
Therefore, these states may be omitted in the evaluation 
of the matrix element of V ra d . Bose statistics furthermore 
require that J' — 1 . Consequently, because of dipole selec- 
tion rules, only the M' = 1 state can be excited starting 
from the M — 2 lower state, which implies cr~ transitions 
induced by laser 1. 

The matrix element involves a spatial integration, which 
can be split up in an angular part and a radial part (the 
so-called Franck- Condon integral). It should be noted that 
the Franck-Condon principle (i.e. approximating the inte- 
gral by replacing it with the integrand evaluated at the 
Condon radius) does not apply to the transition to the 
v = excited state in 0+. This is due to the fact that 
the v = radial wavefunction, unlike more highly excited 
vibrational wavefunctions, has a broad amplitude max- 
imum which spans the entire classically allowed region. 
Also, a common approximation in which only the radial 
wavefunctions u(R) in the FC integral are included is in- 
adequate, since the angular part of the matrix element 



noted that the atomic dipole moment in Eq. J3J) is modi- 
fied so as to take the presence of the magnetic field (which 
lifts the degeneracy of the lower state) into account. 

For the matrix element describing the radiative cou- 
pling in the OFR, the radial integration is limited to the 
internuclear range for which the excited-state wavefunc- 
tion is known ( Section 13.2(1 . thereby neglecting the short- 
range contribution. Since in this case the excited state is 
virtually a pure Hund's case (a) state, the radial integral 
only contains the radial wavefunctions of the lower and 
excited states. 

The bound-bound matrix element in Eq. (TJJ) is ob- 
tained along the same lines as given above, but now in- 
volving an M' = 1 — > M = 2 transition. It is furthermore 
assumed that the magnetic moment of the v = 14 bound 
state in the 5 i7+ potential is identical to that of the initial 
free-atom pair, such that all Zeeman shifts due to the field 
of the magnetic trap cancel out. 
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